Accounting for spin fluctuations beyond LSDA in the density functional theory 
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We present a method to correct the magnetic properties of itinerant systems in local spin density 
approximation (LSDA) and we apply it to the ferromagnetic-paramagnetic transition under pressure 
in a typical itinerant system, N13AI. We obtain a scaling of the critical fluctuations as a function 
of pressure equivalent to the one obtained within Moryia's theory. Moreover we show that in this 
material the role of the bandstructure is crucial in driving the transition. Finally we calculate the 
magnetic moment as a function of pressure, and find that it gives a scaling of the Curie temper- 
ature that is in good agreement with the experiment. The method can be easily extended to the 
antiferromagnetic case and applied, for instance, to the Fe-pnictides in order to correct the LSDA 
magnetic moment. 
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Density functional theory (DFT), in its most com- 
mon implementations (local spin density approximation, 
LSDA, with or without gradient corrections, GGA), is 
in principle the only way to access the ground state of 
real materials X And indeed the agreement with exper- 
iment concerning the ground state properties, such as 
crystal and electronic structures, is excellent, especially 
for itinerant system, where local correlations play a mi- 
nor role. Nevertheless, a well known problem of LSDA is 
the overestimation of the tendency to magnetism in itin- 
erant magnets near the quantum critical point (QCP). 
This problem can be traced down to the fact that LSDA 
is essentially a mean field theory, which does not take into 
account a detrimental effect of near-critical fluctuations 
on the long-range magnetism. Thus, while the itiner- 
ant nature of systems like FeAl, 2 Pd 3 or the more recent 
and better known Fe-pnictides^— make LDA and GGA 
reproduce very well the paramagnetic bandstructure, 
whenever a (magnetic) quantum critical point (QCP) is 
approached the theory fails miserably. The importance 
of this problem is demonstrated by the amount of papers 
dealing with the problem of correcting the magnetic mo- 
ment of Fe-pnicitides^r— There, the usual argument is 
that correlations beyond mean field suppress the (LSDA) 
local ordered moment. It was shown that one can reduce 
the calculated magnetic moment by using the LDA+U 
method with a negative U», but there is no physical 
justification for this procedure. Such many-body ap- 
proaches as Dynamic Mean Field Theory (DMFT)2r— 
and Gutzwiller— were also successfully used; since these 
methods introduce additional fluctuations into the sys- 
tem, they obviously work in the right direction. How- 
ever, the concept of substituting long-range critical fluc- 
tuations by the on-site ones is rather questionable^ Fur- 
thermore the effect of non-local fluctuations was recently 
found to be crucial, also in localized models, whenever 
the critical behavior is analyzed—, and in any event com- 
putational load in these methods is incomparably heavier 
than in LDA calculations. 



For these reasons we propose a different approach 
which corrects LSDA within DFT and takes into ac- 
count the itinerant nature of the system. Our method 
is easy to implement^ carries no additional computa- 
tional cost, and has a transparent physical justification. 
The approach we describe in the following is based on 
the idea that unknown more accurate DFT is very close 
to the conventional LSDA-GGA functional, but the en- 
ergy gain due to spin polarization ( "Stoner interaction" , 
using the DFT parlance) is reduced by about as much 
as the Moriya selfconsistent renormalization (SCR) the- 
ory^, successfully used befor o 3 ' 16 ' 17 , suggests. For this 
reason we call it reduced Stoner theory (RST). In fact we 
show that there is a one to one and well defined connec- 
tion between our method and the Moriya SCR theory in 
accounting for the effect of spin fluctuations in itinerant 
magnets. This allows us to make an ab-initio predic- 
tion of the magnetic moment as a function of pressure 
for the archetypical N13AI itinerant ferromagnet, with 
the correct scaling of the Curie temperature, which un- 
til now was impossible. For M3AI indeed, the effect of 
spin fluctuations on physical properties both in the mag- 
netic state (magnetic moment and susceptibility) and in 
the paramagnetic one (dc resistivity) was demonstrated 
both theoreticall y 16 ' 18 and experimentally! 19 ' 20 

Applying the RST we estimate the pressure depen- 
dence of the contribution due to spin fluctuations on the 
magnetic moment of N13 Al. Surprisingly we find that this 
contribution is almost pressure independent. This means 
that spin fluctuations only act in shifting the Stoner con- 
dition. The rest is done by the bandstructure. As a 
consequence we show that the way in which the spin 
fluctuations renormalize the effective Stoner interaction 
is encoded in the bandstructure itself. This is generally 
true but in M3AI it becomes particularly evident due to 
the perfect scaling of the density of states (DOS) with 
pressure. 

Having established a link between the RST and the 
SCR theory indeed, we are able to give a reliable estimate 
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FIG. 1: Top: Magnetic moment of M3AI per Ni atom calculated 
in LSDA (black dotted line). Both the magnetic moment at zero 
pressure m(0) and the critical pressure P c are overestimated with 
respect to the experimental ones (red triangles). The green dots 
mark a linear interpolation of the data. Bottom: Paramagnetic 
density of states as a function of energy calculated for different 
pressures in LDA. Inset: Density of states at the Fermi level calcu- 
lated as a function of pressure (black squares). The data are fitted 
with eq. [TJ(red continuous line). 



of the Stoner parameter in itinerant magnets. Moreover 
if one reverses the logic a comparison between the ex- 
periment and the RST results gives an easy and reliable 
estimate of the spin fluctuations acting in the system. 
In the following we present the bare LSDA results for 
N13AI under pressure and we interpret them within the 
so-called extended Stoner Theory (EST) We explain 
that the overestimation of both the magnetic moment at 
zero pressure, and the critical pressure P c come essen- 
tially by the overestimation of the Stoner parameter in 
LSDA. After that we introduce the formalism for cor- 
recting the LSDA behavior and we present the scaling 
equations used in the RST. In the last part of the paper 
we apply these scaling equations to the ferromagnetic- 
paramagnetic transition of M3AI under pressure. 

First we discuss the bare LDA (and LSDA) results for 
the magnetic transition in M3AI under pressure. M3AI 
crystallizes in the ideal cubic CU3AU cP4 structure. We 
calculated the equilibrium lattice parameter a = 3.4825 
A in LDA 24 and we found it to be ~ 2% smaller than the 
experimental one^ The calculated magnetic moment at 
P = GPa is m(0) = 0.68 fi B in LSDA (0.73 fi B in 
GGA) in reasonable agreement with previous results^. 
After that we calculated the magnetic moment as a func- 
tion of pressure and, as shown in the top of Fig. [TJ we 
found that it decreases approximately linearly up to the 
critical pressure P c = 45 GPa, much larger than the ex- 
perimental one, which is P c Ex p ps 8 GPa,- 20 This behavior 
can be easily understood within the EST* 2 ^— In fact, as 
shown in Fig. [1] (bottom), the paramagnetic DOS N(E) 



calculated in LDA scales almost perfectly with the pres- 
sure P as: 



N(E,P) 



N(E*Z(P),0) 

W) : 



(1) 



in a wide energy range around the Fermi level taken as 
the origin, where Z(P) = 1 + 0.005 * P with P measured 
in GPa. This can be related to magnetization via the 
EST— ~— , which combines the Stoner criterion with the 
Andersen's force theorem to show that in the lowest order 
in the magnetic moment m the total LSDA energy is 
given in terms of the non magnetic DOS per spin: 



E(m) 



m! dm' 1 

— Ira 1 

N(m>) 4 



(2) 



where iV(m') is the extended DOS defined in Ref. |2l| as 
the average DOS over an energy interval equal to the ex- 
change splitting at a given m, and I is the Stoner param- 
eter. Stationary solutions appear where N(m) = 1/7, 
and dN{m)/dm < 0, forming a stable or a metastable 
(metamagnetic) ferromagnetic state with the magnetic 
moment m. If N(E) scales according to eq. [TJ so does 
also N(m). Comparing eq. [5]with the fixed spin moment 
calculation we found the Stoner parameter in LSDA to 
be I LSDA ss 0.41 eV independent on pressure (see ap- 
pendix |XJ). In GGA the Stoner parameter is 17% larger, 
jGGA _ i.i-jjLSDA j n tne f ii ow i n g we w m consider 

only the LSDA results. Around the minimum of E(m) 
N{m) decreases monotonically. Thus a reduction of 
I shifts the Stoner condition [N(m) = 1/1] toward a 
smaller magnetic moment, bringing the LSDA results in 
agreement with the experiment. The aim of this paper 
is to find a way for reducing the LSDA Stoner parameter 
and getting in this way the correct magnetic moment in 
itinerant systems. In order to gain some understanding 
of how the standard LSDA needs to be corrected to ac- 
count for spin fluctuations, let us recall the SCR theory 
for critical ferromagnets. It starts with the Ginzburg- 
Landau expansion, 



E(M) =a Q + a 2 M 2 /2 + ai M 4 /4: + a 6 M 6 /6 + 



(3) 



where M is the magnetic moment and E is the LSDA 
energy. It is then assumed that M fluctuates around the 
average value M, so that M = M + SM. Assuming that 
SM follows a Gaussian distribution such that = 
£ 2 , we can rewrite eq. [3] in terms of M as 



E(M) = a + & 2 M 2 /2 + d 4 M 4 /4 + d 6 M 6 /6 



(4) 



where the explicit formula for the renormalized coeffi- 
cients is given in Ref. y, and the second order coefficient, 
a 2 = a 2 + (5/3)a 4 £ 2 + (35/9)a 6 £ 4 + • • • . If we restrict 
our expansion by the second order, this SCR procedure 
is equivalent to renormalizing spin susceptibility. 

The amplitude of the spin fluctuations, £, in principle, 
can be obtained by the fluctuation-dissipation theorem 
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(FDT): 



set of scaled equations: 
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f/q 
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7mx(q,w) 



(5) 



where f] is the volume of the unit celli^ xili w ) m this 
formula is often expressed in terms of the noninteracting 
(Lindhard) susceptibility, which in turn is expanded to 
the lowest order in q and u), xo(q, w) = xo(0, 0) — aq 2 + 
ibuj/qZ. The coefficients a and b can be written down as 
functions of the Fermi velocity, averaged over the Fermi 
surface. However, in order to calculate £ using this ex- 
pression one needs to apply a cutoff q c which is not a well 
defined quantity. 

Let us now point out that the uniform spin suscep- 
tibility in the LSDA can be written in a particularly 
simple form, namely x _1 = S 2 E/6M 2 |m=o= a i — 
\ n (E F ) ~ I (this can be considered as a rigorous def- 
inition of the LSDA Stoner parameter I), where N a {Ep) 
is the paramagnetic DOS per spin. Comparing this ex- 
pression with the one for 02 we see that the SCR proce- 
dure is equivalent to renormalization of / according to 



1 = 1— (10/3)a 4 £ 2 



(6) 



which can also be written as I = si, where s < 1 (see 
also appendix |A"| . 

This is a justification of the recipe of using / as an ad- 
justable parameter—, often used in critical ferromagnets 
empirically In this sense, / can be perceived as derivable 
from an unknown more accurate DFT, in a specific ma- 
terial. One can also "reverse-engineer" such an improved 
functional, using the standard von Bart-Hedin scaling, 



E r 



e xc (n,m)n(r)dr 



£xc(n,0 = £xc( n ) + /(C)Ae a: c(n), 



(7) 
(8) 



where £(r) = m(r)/n(r), e xc and Ae xc (n) do not depend 
on m, and /(C) is a known function, and n = (n-f + n^), 
m = (n-f — ni). The response to magnetism is entirely de- 
fined by the Ae xc (n) functional, as the energy difference 
between the fully polarized and unpolarized electron gas: 



de xc 
dn 
de xc 



del 



dn 



on n 



= /'(OAMn)- 
am n 



(10) 



Note that the charge potential also acquires a term that 
disappears when £ = 0. In this sense, it is impossible 
in a physically meaningful way to reverse-engineer / and 
Ae xc in such a way that de xc /dm be scaled by a constant 
factor s, and de xc /dn would remain the same. Rather, a 
natural way to weaken the magnetism in this formalism is 
to scale Ae xc {n) in eq. |S] Then we will have the following 



Sxc(n, C) = £ xc (n) + sf(QAs xc (n) 



(11) 



= < c + s[f(C)Av xc (n) - f'(0As xc (n)-] (12) 
on n 



om n 



(13) 



where the part in the square brackets is simply the ad- 
ditional charge potential that appears because of spin 
polarization. It is easy to verify that this functional pro- 
duces an exchange-correlation potential scaled by s, and 
the charge potential unchanged: 



Vt(r) 



-n(r) 



S [F t (r)-n(r)] 
[F t (r)+^(r)] 



(14) 



and the Stoner kernel 5 2 E XC /Sm(r) 2 also scaled by s, 
as we wanted. Eq. [14] can be easily implemented^ and 
can be used to obtain correct magnetic moments and the 
corresponding electronic structure in the materials near 
fcrro or anii/erromagnetic QCP. Moreover, given eq. [5] 
and [6] s gives also an indication of the strength of spin 
fluctuations acting in the system. Below we use eq. [TJ] 
for correcting the LSDA results obtained in the previous 
section in the case of M3AI itinerant ferromagnet under 
pressure. As shown previously, N13AI becomes paramag- 
netic under pressure^ and LSDA overestimates not only 
the magnetic moment, but also the critical pressure P c . 
Empirically, by comparing the LSDA value of the mag- 
netic moment with the experimental one for P = GPa 
and for P = 6 GPa, we found that in both cases the value 
of s needed to reconcile the LSDA result with experi- 
ment, using the scaling introduced in eq. HH is s « 0.88. 
This implies that £ is almost pressure-independent be- 
tween and 6 GPa. This value of s gives for I the same 
value obtained by renormalizing / within the EST (see 
appendix [A]) . In general, one expects that spin fluctua- 
tions become stronger closer to the critical pressure, so 
that their average amplitude £ becomes larger, and the 
scaling parameter s = 1 — (10/3)a4£ 2 // in eq. [5] smaller. 
On the other hand, £ is defined by averaging over the en- 
tire Brillouin zone, and the fact that susceptibility at one 
particular point q = diverges may or may not strongly 
affect £. In order to understand this result, we compared 
our calculations with Ref. [T7J, where £ is calculated for 
P = GPa in the approximations described above. By 
means of scaling arguments described below, we found 
that the transition is driven entirely by the change in 
DOS given by eq. Q] while the renormalized Stoner pa- 
rameter / is, in the first approximation, pressure inde- 
pendent. Let us try to rationalize this fact, using the 
fluctuation-dissipation theorem (eq. [5])- 

Given eq. [1] and since the Fermi velocity scales in- 
versely with the DOS, v(P) = v(0)Z(P), in the approxi- 
mation described in the beginning, and used in Refs. [ill 
we found that £ does actually scale with pressure as 



£(P) a 



TttP) 



where Q(P) is the unit cell volume at 
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analyzed the data using a Landau functional for the field 
H = a 2 m + 04m 3 - cV 2 m. (15) 
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FIG. 2: Predicted magnetic moment as a function of pressure 
calculated within the RST (red dotted line). The critical behavior 
in oc (P c — P) 1 / 2 (blue line) is followed between 2 and 6 GPa. 
The Curie temperature T c oc m 3 / 2 is shown in the inset by black 
dots. Apart from the underestimation of the critical pressure, the 
agreement with the experimental data taken from Ref. (green 
dots) is very good. In fact, both theoretically and experimentally, 
the extrapolation of the data the data at 5.7 GPa (black dashed 
line) and 8 GPa (green dashed line) respectively, give both the 
same critical behavior T c oc (P — P c ) 3 / 4 , as observed by Nicklovic 
et al22. 



pressure Pi^L We found indeed that it does not diverge or 
even grow substantially near the critical pressure (even 
though the model susceptibility at q = 0, u = does 
diverge there, IN(Ep) — 1). This result is particu- 
larly important because it tells us that the effect of spin 
fluctuations in N13AI under pressure can be entirely ac- 
counted for by renormalizing the Stoner factor (the rest 
is a bandstructure effect). Furthermore, the fact that 
the correlation length associated with the fluctuations (£) 
does not diverge approaching the transition, is compat- 
ible with the fact that the system does not show other 
instabilities (like triplet superconductivity) at P = P c . 
In fact if £ would go to infinity, another (competing) in- 
stability could profit of this kind of long range correla- 
tion in order to build up a competing order parameter. 
In fact, yjn(6 GPa)/ft(0 GPa) = 0.98, which implies 
that £ changes only by 2%. Applying the SCR theory 
starting from the fixed spin moment calculations and us- 
ing £ as a parameter (see appendix [A} , we found that 
the best agreement with the experiment was achieved 
if £ changes slightly more, by 5%, but simply using a 
pressure-independent £, corresponding to the scaling pa- 
rameter s — 0.88, provides, apart from some underes- 
timation of P c , a very reasonable agreement with the 
experiment. This choice of s allows us also to make pre- 
dictions about the magnetic moment between and P c , 
as shown in Fig. [5J 

There are no experimental measurements of magneti- 
zation as a function of pressure between and 8 GPa in 
N13AI, but the pressure dependence of the Curie temper- 
ature T c was measured by Nicklovic et al£2- They also 



By assuming a 2 oc (P — P c ) they found T c oc (P — P c ) 3//4 
in good agreement with the experimental data. Within 
Moriya's theory^ T c oc m 3 / 2 where m — \ m |. Assuming 
a 2 oc (P-P c ) in eci.lTSl gives m oc (P-P c ) with ^ = 1/2, 
from which the result of Nicklovic et al. follows. In the 
following we show that the linearity of a 2 with respect to 
(P — P c ) is a consequence of eq. [TJ For small values of the 
magnetic moment, in the fluctuation-corrected LSDA de- 
scribed above, a 2 = 2N \e f ) — i> c = ^> an< ^ ? ^ s adjusted 

1 at P = P c . Given eq. Q]thcn 
' [a(P-P c )\ cx (P-Pc), 



so as to have N a (Ep)I 

_ Z(P)/Z{P C ) i _ 

u>2 



2N„(E F ),P C ) 2 _ 2Z(P C ) 

thus providing a microscopical justification for the model 
of Ref. 20. As shown in Fig. [T] in LSDA (for large values 
of the magnetic moment) we found j3 w 1, while j3 = 0.5 
is recovered for small values (see Fig [5J . This is due to 
the fact that at large magnetic moment the coefficient a 2 
must be corrected by adding high order terms. Finally, 
using Moriya's relation for T c we find T c oc (P — Pc) 3 / 4 
in full agreement with Nicklovic et al. The disagreement 
with experiment in Fig. [2] (inset) concerns only the un- 
derestimation of P c caused by: the approximation of s as 
a constant value and the underestimation of the equilib- 
rium lattice parameter^ 

To summarize, in this paper we have described a 
simplified method for accounting for near-critical spin- 
fluctuations within the DFT. The method amounts to 
scaling the DFT exchange-correlation field by a phe- 
nomenological constant s, and subsequent self-consistent 
solution of the Kohn-Sham equations. This phenomeno- 
logical constant can also be, in principle, calculated via 
the fluctuation dissipation theorem 3 -, and in this sense is 
equivalent to the SCR theory by Moriya. Our method 
is complementary to the widely used similarly semi- 
phenomenological LDA-t-U, and plays for itinerant near- 
critical magnets the same role as LDA+U for systems 
near a Mott-Hubbard transition. We apply this method 
to the ferromagnetic QCP in M3AI under pressure. We 
show that, due to a particular scaling property of the 
bandstructure, the parameter s is constant with pres- 
sure. In this way the method becomes completely ab 
initio. In fact, the ferromagnetic-paramagetic transition 
in NisAl is driven by the band structure changes under 
pressure, while the feedback to the critical fluctuations 
(parameterized by s) is small. Using this formalism, we 
make a prediction of the magnetic moment as a function 
of pressure, which recovers the critical exponent for the 
magnetization j3 = 1/2 for small magnetic moments and 
explains the observed dependence of the Curie tempera- 
ture on pressure. The new method should be useful in 
cases when one needs to calculate electronic properties 
of materials where LSDA overestimates the tendency to 
magnetism, or when one wants to monitor theoretically 
the evolution of the electronic structure from non mag- 
netic to magnetic and get an estimate of the spin fluctu- 
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Appendix A: Connection between RST and SRC via 

EST 

If we fit the fixed spin moment calculations with a Lan- 
dau expansion up to the 8-th order we find the same re- 
sult as in Ref. 161 However, as opposed to that work, 
we do not want to apply explicitly a fluctuation-induced 
renormalization, but renormalize / directly, according to 
eq. 1141 that is, to scale down ab-initio the Stoner param- 
eter in order to obtain the correct magnetic moment for 
N13AI. To make a direct connection between the Landau 
approach and the RST approach, we can use the EST. 
Indeed, as discussed above, the RST method is equiva- 
lent to the SCR method, if only the second order coeffi- 
cient ai (the inverse magnetic susceptibility) is renormal- 
ized. We can expand the EST total energy in eq. [5]in to, 
and then renormalizing the Stoner parameter becomes 
exactly equivalent to a renormalization of a?.. Fig. [3] 
shows that, in fact, only the ai renormalization affects 
the equilibrium magnetic moment, while the renormal- 
ization of the higher order coefficients influences E(m) at 
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FIG. 3: Energy E as a function of the magnetic moment m for 
P = GPa. The red dots show the curve obtained by means of 
EST with a renormalizcd Stoner parameter I = 0.358. The green 
and the blue lines mark respectively the Landau functional where 
all the coefficients are renormalized and the one where only the 
first coefficient of the expansion is renormalized. 



m larger than the experimental magnetic moment. Note 
that the red dotted curve (EST) in Fig. [3] corresponds to 
rigid splitting of the unpolarized bands^i, while the coef- 
ficients of the Landau expansion were obtained by fitting 
fully selfconsistent magnetic calculations. 1 ^ In Fig. [4j is 
shown the equilibrium magnetic moment m as a function 
of the reduced Stoner parameter / = si, obtained both 
in the EST and by self-consistent calculations using the 
exchange scaling s of eq. UH at P = GPa. The value 
of I is chosen, in the former case, in such a way that the 
LSDA and EST magnetic moments have the same value 
for s = 1. The agreement between the two curves is per- 
fect and the value of I found in this way is in perfect 
agreement with the one extracted from the comparison 
of the EST with the renormalized Landau expansion in 
Fig. El 
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FIG. 4: Magnetic moment m as a function of the reduced Stoner 
parameter I = s * I for NiaAl at GPa in EST (red dots) and 
RST (green dotted line). The bare value of te Stoner parameter 
I = 0.407 eV was chosen in order to to have the same value of m 
for s = 1. 



